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Abstract 

IceCube has measured a diffuse astrophysical flux of TeV-PeV neutrinos. The 
most plausible sources are unique high energy cosmic ray accelerators like hypernova 
remnants (HNRs) and remnants from gamma ray bursts in star-burst galaxies, which 
can produce primary cosmic rays with the required energies and abundance. In this 
case, however, ordinary supernova remnants (SNRs), which are far more abundant than 
HNRs, produce a comparable or larger neutrino flux in the ranges up to 100-150 TeV 
energies, implying a spectral break in the IceCube signal around these energies. The 
SNRs contribution in the diffuse flux up to these hundred TeV energies provides a 
natural baseline and then constrains the expected PeV flux. 



1 Introduction 


The IceCube (IC) Collaboration has found evidence of a flux of extraterrestrial neutrinos 
in the analysis of the 3-year combined data [1, 2 , 3]. Although, the exact origin of these 
events is yet to be settled, the signal seems to be incompatible with the atmospheric neutrino 
background and only a small fraction of these events is consistent with a galactic origin [4] . 
The flux can be reasonably considered to be nearly isotropic [2] and can be interpreted as 
the ‘smoking gun’ signal for hadronic processes [5, 6 ] in the cosmic ray (CR) accelerators. 
Indeed, such neutrinos would also put extremely stringent constraint on neutrino Lorentz- 
invariance violation [7]. 

Several studies connecting these IC neutrino events with different cosmic ray produc¬ 
tion sources estimated the diffuse TeV-PeV background (see [ 8 ] for review and references 
therein). The general idea being that the CRs loose energy via the hadronic processes (py 
or pp collisions) with the resultant meson decay to generate neutrinos and gamma rays. 
However, the explanations are based on several assumptions and free parameters, also due 
to the poor statistics of events too many explanations seem compatible and is not leading 
to any reasonable conclusion. 

Among these possibilities the explanations that these events are connected to the diffuse 
background of neutrinos from extreme high energy accelerators like semi-relativistic hyper¬ 
nova remnants (hereafter ‘HNRs’) [9, 10, 11], long GRB remnant [12] embedded in giant 
molecular cloud or slow jet active galactic nuclei (AGN) [13] with dominant contribution 
from star forming galaxies seem to be highly plausible. Such objects in star forming galaxies 
can be tuned to explain the IG diffuse neutrino flux normalization and the characteristic 
sharp cut off at PeV energies. In particular, the 7 rays generated in the same hadronuclear 
processes would also populate the diffuse gamma ray background and the observed Fermi 
diffuse background [14] can constrain these models [11, 13]. 

Among the stellar remnants supernova remnants (SNRs) can accelerate protons to PeV 
energies [15], whereas the HNRs and GRB remnants, due to their greater explosion energy 
and shock velocity, can generate protons up to hundreds of PeV [16] . The secondary neutrinos 
carry almost 5% of the parent proton energy, implying that the normal SNRs can contribute 
to neutrino fluxes in the 100 TeV energy range whereas the more energetic GR accelerators 
may generate neutrinos to PeV energies. The far away, metal poor galaxies with large star 
formation rate and large proton densities are considered to be strong sources of such remnants 
[17]. In addition, the large pion production efficiency in these star burst galaxies (SBGs) 
makes them excellent sources of the diffuse secondary neutrinos [18]. The rate of the SBGs 
increases with the redshift and will have the dominant contribution from high redshifts {z ~ 
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1-2). Therefore one needs to consider the redshift integrated diffuse neutrino spectra, which 
depends on the proton accelerating power of the remnant in the host galaxy and the rate of 
such remnants. The normal SNRs are expected to contribute to this diffuse neutrino flux up 
to the 100 TeV energy range and the remnants with greater ejection energies in the whole 
TeV-PeV range. 

However these extreme high energy remnants like hypernovae and GRBs are rare events, 
with a local rate of less than 1% of the normal SNRs [19, 20]. The SNRs being 100 times 
more abundant would have larger or at least similar energy budget than the HNRs. Thus the 
normal SNR component of the diffuse neutrinos in the 100 TeV energies should be at least of 
the same order as the HNRs. Also the fact that IC is observing reasonable number of events 
below the 100 TeV energies [3] emphasizes the importance of the normal SNR contribution 
to the diffuse flux. 

In the following, we calculate the diffuse TeV neutrino flux for SNRs. For the TeV-PeV 
diffuse flux we take the example of the HNRs. Our calculations show that for a similar popu¬ 
lation of host galaxies the neutrino flux up to 100-150 TeV is dominated by the so/fer spectra 
of the SNRs component and also has a substantial contribution from the HNRs. Whereas, 
the higher energy part (100-1000 TeV) is only HNR generated and has a harder spectra, 
with a steep cut off around PeV energies. Therefore, the diffuse flux would have different 
spectral behavior in different energy ranges, implying a spectral break around 100 TeV. The 
present statistics of the IC data is insufficient to conhrm or exclude this picture and future 
IC data should test this energy dependent spectral behavior. Also, the dominant flux at 
the 100 TeV energies will give the normalization of the diffuse flux and can put interesting 
constraints on the connection between SNR and HNR rate. 

The paper is outlined as follows. In section 2 we describe the diffuse neutrino background 
from SNR and HNR. We give a brief overview of the uncertainties which go into the calcu¬ 
lation. In the next section 3 we show the results and compare the flux with the IC events. 
Finally in section 4 we conclude with a review of our results. 


2 Diffuse neutrino background from stellar remnants 

The observed diffuse flux of neutrinos from a particular stellar remnant would have contri¬ 
butions from different redshifts, 

dN{E:^) ^ p- dN{E,) Rsr{z) dz 

dEf AttHo Jo dE, + ^)='+ Ca ’ 
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where the -Rsr(^) is the rate of the stellar remnant (SNR or HNR) and dN[Ey] / dE^ is 
the source neutrino flux. The observed energy {E'^) is connected to the original energy 
{Ey = {1 + z)Ey^). The Hubble parameter (Hq) used is 0.69 km s“^Mpc“^ For the standard 
A-CDM cosmology, the matter and dark energy density Hm and Ha are taken to be 0.27 
and 0.73, respectively [21]. The source flux dN(Ey) / dEy depends on the host galaxy and 
thus would depend on the relative rate of the different types of galaxies. In particular, the 
relative population of the SBGs (/sbg) is estimated to be 0.1-0.2 [22, 23, 24] and the rest 
is considered to be normal star forming (NSFG) ones. Thus, total flux from a particular 
remnant population is the weighted sum of flux from SBGs and NSFGs. 


The contribution form each stellar remnant from a particular type of host galaxy is given 
by [25], 


dN{Ey 
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where is efficiency of the pion production, k is the inelasticity (0.2), Ey is the secondary 
neutrino spectrum and Jp is the primary proton spectrum Ep^exp{-Ep/E^^^). The 
normalization of the proton spectrum is estimated from the total proton energy {Ej), which 
is a fraction of the ejected energy of the stellar remnant. In particular, for the HNRs Ej is 
in the range 5 x lO^^-lO^^ erg [26], whereas for SNRs the total proton energy is expected to 
be at least one order lower. The maximum energy (E™®'^) up to which the parent protons 
can be accelerated is governed by the total ejecta energy, magnetic field and the shock radius 
of the remnant. The HNRs with larger ejecta energy and shock radius are expected to have 
stronger neutrino production with broader spectra reaching much higher energies compared 
to the normal SNRs. The for SNRs and HNRs are considered to be in the range 1-10 
PeV and 10^-10^ PeV, respectively [15, 16]. 


The pion production efficiency (r^jr) depends on the properties of the host galaxy envi¬ 
ronment [11]. Efficient pion production would require that the energy loss time (tioss) due to 
pp collisions is smaller than the proton conhnement time (tconf), thus the is estimated as 
? 7 ^ = min(l, tconf/tioss)- Whether tconf is long enough (> tioss) to loose energy via collisions 
with the interstellar medium (ISM) gas would depend on the gas density and the magnetic 
field strength of the galactic environment. In particular, tioss = [KO'ppUpc]”^, where Upp is 
the inelastic nuclear collision cross section [25]. Glearly, because of this inverse dependence 
on the (up), larger proton dense environment will have smaller energy loss time. There¬ 
fore, SBGs, with larger proton density (up few lO^cm ^) [27] compared to the NSFGs 
(up ~ 10 cm“^) [28, 29], are expected to have a more efficient energy loss. 

The tconf depends on the scale height of the galaxies (h), the diffusion coefficient (D) and 
the galactic wind velocities (Vwind)- At low energies tconf is dominated by the advective 
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escape [30] via galactic wind (tconf = tadv = h/Vwind) and in the higher energies diffusive 
escape takes the lead (tconf = tdifr = h^/4D) [31]. The energy range (Eb) where both 
processes become competing would result in a break of the spectra and for energies above 
this break the diffusive escape would give a softer spectra with a cutoff-like feature. In 
particular, the effective spectral shape of the neutrinos in the advection dominating regime 
broadly follow the primary proton spectra At higher energies the energy dependence 

of the diffusion coefficient {D oc E~^'^) would make the resultant neutrino spectra softer. In 
addition, at high energies, the app would also have an effect on the overall neutrino spectral 
shape. 

The break energy (Eb) would also depend on the galactic properties. For SBGs the break 
would appear at much higher energies (PeV) compared to the NSFGs (few hundred GeV) 
[10] implying the fact that the fluxes in the TeV-PeV energy range would be dominated by 
SBGs. In the advective escape regime the galactic wind velocities in SBGs are much larger 
than the NSFGs due to the higher rate of supernova explosions, so that one may expect a 
value of 1500 km s“^ and 500 km s“^ for SBGs and NSFGs, respectively. Regarding the 
diffusive escape, the high redshift SBGs are expected to have a lower diffusion normalization 
(10^^cm^s“^) because of the higher magnetic helds, whereas for the NSFGs we assume a 
similar diffusion normalization (10^®cm^s“^) to the one of our galaxy [32]. Also the scale 
height of the NSFGs (~ Ikpc) is larger compared to the SBGs (~ 0.5kpc). Adding all these 
facts, the overall pion production efficiency of the protons is much larger in SBGs, resulting 
in larger neutrino fluxes compared to the NSFGs. For example, at 100 TeV the is around 
0.02 and 0.6 for the NSFGs and SBGs, respectively. 

The total high energy diffuse neutrino background would also depend on the relative pop¬ 
ulation of the different stellar remnants. In particular, the Rsr{z) follow the star formation 
history i?sFR(2;), for normal SNRs RsNuiz) = 1.22 x 10“^i?sFR(2:)MQ^ [33]. The local hyper¬ 
novae rate is about 1% of the supernovae rate [19, 20, 34] and in our analysis we consider 
a HNR rate i?HNR(^) < 10 “'^Rsfr(2;)Mq^ For the Rsfr(^) we use the concordance models 
of [35, 36], where in low z (< 1) the SFR increases as (1 -|- to remain constant till the 
redshifts x ~ 4. Beyond z = 4 the SFR decreases in extreme rapid fashion ~ (1 -|- and 
for the local star formation rate i?sFR(0), we use 10“^MQyr“^Mpc“^. 


3 Results 

The flux of high energy diffuse neutrino background depends on the considerations described 
in the previous section. To study the dependence on the uncertainties of the input parameters 


5 



Figure 1: Single flavor diffuse neutrino spectra from both the SNRs (red continuous curve) 
and HNRs (blue dashed) contributions. The estimates have contribution form both the 
normal and the star burst galaxies. The band describes the uncertainties in the input 
parameters. The IC best £t spectra and the data is given by the black continuous curve and 
black points, respectively [2]. 

we calculate the flux for the extreme ranges of the parameters. In particular, the input 
parameters like total proton energy (ifj) and the relative population (/sbg) of the SBGs in 
comparison to the NSFGs are difficult to estimate accurately. For the HNRs we use E'^ and 
/sbg in the range 5 x erg and 0.1-0.2, respectively. In case of the SNRs the 

is considered in the 5 x erg range. Also, there will be uncertainties in the local 

HNR rates. The value used in the analysis is simply the upper limit. However, in the SBGs 
the HNRs may have a higher rate compared to the normal star forming galaxies [36, 37]. 
Regarding the we take 5 PeV and 1 EeV for the SNRs and HNRs, respectively. 

Fig. 1 shows the single flavor diffuse neutrino spectra {E^dN/dE^) for both the SNRs 
(red continuous curve) and HNRs (blue dashed) contributions. For both the SNR and HNR 
curve /sbg is taken as 0.15. The Ej is 5 x 10®° erg and 5 x 10®^ erg for the SNR and 
HNR, respectively. The bands around the curves describe the uncertainties in /sbg and Ej. 
The IG data with error bar is given by black points and the best fit spectra {EldN/dEy = 
1.5 X 10“^(E/100TeV)“°'° GeVcm“^s“^sr“^) [2] is given by the black thick curve. Both the 
HNR and SNR estimates have contributions from the normal and the star burst galaxies. 
Due to high and Eb in the TeV-PeV energies the SBGs contribution dominates over the 
normal ones. The HNR contribution steeply decreases around few hundred TeV to give a 
sharp cut off around PeV energies. This feature of the diffuse HNR neutrino flux makes it 
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Figure 2: The total diffuse neutrino ffux from both stellar remnants is given by the purple 
curve. The band describes the uncertainties in the input parameters. The IC data is the 
same as in fig 1. 

an excellent candidate to the observed IC events extended to PeV energies. However, in the 
lower TeV energies the same galaxy population is also giving a strong diffuse flux coming 
from the normal SNRs. The SNRs contribution is dominating up to hundred TeV energies. 
The limitation of the normal SNRs to accelerate protons above PeV energies translates into 
this cutoff at around 100 TeV in the secondary neutrino spectra. 

The total diffuse neutrino flux from both types of stellar remnants and all galaxy pop¬ 
ulation is given in the Fig. 2. The purple curve is the total diffuse flux corresponding to 
the SNR and HNR curves in Fig. 1. The band around the total flux curve describes the 
uncertainties in the input parameters. The IC best fit and the data are again plotted by the 
black curve and the points, respectively. Up to the 100 TeV energies the spectral shape of 
the diffuse neutrino flux is dominated by the SNRs contribution. As the SNRs contributions 
rapidly decrease around 100 TeV, the higher energy contribution is solely from the HNRs. 
The HNRs flux in the higher energies is harder compared to the softer SNR contribution. 
Thus the total diffuse neutrino spectra will have a break in around 100 TeV energies. This 
structure in the spectra would be a generic feature of these models and can put interesting 
constraints. Indeed, the uncertainties in the model parameters can only effect the overall 
normalization of the diffuse flux, but the spectral structure remains similar . 

In particular, the extreme high energy accelerators like HNRs have many uncertainties. 
The shock acceleration, the magnetic fields and efficiency of kinetic energy to CR energy 
transfer in these rare objects still need detailed understanding. The relative rate of HNR 
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to SNR is another ingredient in these nncertainties. However, the SNRs are much better 
understood and their rate is also well estimated. In fact, the MeV neutrinos emitted during 
the supernova explosion phase of these objects would contribute to the diffuse SN neutrino 
background (DSNB). Together with the peculiar neutrino oscillations in SN, DSNB is one 
of the very interesting signals to study [39, 40]. The Super Kamiokande (SK) upper limit 
of DSNB puts a strong limit on the SN rate [38]. Thus compared to the uncertain HNR 
generated diffuse flux the SNR contribution can be considered as ‘guaranteed’ component of 
the diffuse high energy neutrinos. For any high energy neutrino observation this dominant 
and guaranteed SNR contribution would be the natural baseline to hnd the diffuse flux 
normalization in the sub 100 TeV energies. In fact, as pointed out in Fig. 2, extending 
these normalization to higher energies will put strong limits on the HNR contribution and 
they may become weaker to explain the IC events in the PeV energies. The observed IC 
data shows a mild preference in the low energies. The present IC flux would allow a large 
fraction of the input parameters in the PeV energies as the flux quickly decreases due to 
diffusion dominated CR acceleration. In the low energies up to few 100 TeV the diffuse 
neutrino parameters are constrained by the present IC observed flux. However, the present 
statistically limited IC observation can not conhrm or discard these models. Future IC 
data should reveal if any such spectral feature exists in the isotropic neutrino flux. With 
the present data one can only point out the strong connection between the fluxes of the 
TeV-PeV diffuse neutrino background in different energy ranges. 


4 Conclusions 

Very high energy cosmic rays accelerated by the stellar remnants of objects like semi rela¬ 
tivistic hypernovae and relativistic GRB remnants interacting with giant molecular clouds 
through hadronic processes are expected to produce diffuse TeV-PeV neutrino flux within 
the reach of IC sensitivities. We found that indeed, the detected IC neutrino flux can contain 
such an explanation with diffuse neutrino flux originating in stellar remnants in high redshift, 
metal poor star burst galaxies. However, the normal SNR population in such galaxies also 
contribute to this diffuse neutrino flux. Our estimation of such a diffuse neutrino flux shows 
that up to 100-150 TeV energies, the diffuse background is dominated by the contribution 
from normal SNRs. In higher energies the diffuse flux is dominated by the HNR component. 

Therefore, such a HNR model of diffuse high energy neutrinos successfully explaining 
the IC PeV events would have strong constrain from the SNRs contribution form the TeV 
energies, so that it doesn’t overpopulate the flux below 100 TeV. The present IC data 


indicates mild tension with a large part of the parameter space of these diffuse models. 
Future IC data should give a conclusive indication from the spectral shape of the flux. In 
fact the results would also depend on the connection between the star formation history of 
the different remnants. In the high redshift galaxies the very high energy accelerators are 
expected to be more abundant compared to the local galaxy population, the absence of such 
a spectral shape would constrain the relative rate and redshift dependence of the different 
remnants. Such a detailed study would give more exhaustive answers and should be pursued 
in future. 

Furthermore, the same hadronic interactions of the stellar remnant CRs would also gen¬ 
erate a diffuse gamma ray component from the neutral meson decay and would populate 
the observed diffuse gamma ray background in the 10-100 GeV energies. The correspond¬ 
ing diffuse gamma ray background of the detected IC neutrino flux reasonably agrees with 
the Fermi limits [10, 13]. Thus, such a model of stellar remnants producing the diffuse 
neutrino spectra in star forming galaxies would also get constrained from the Fermi limits. 
The ordinary SNRs are more susceptible to overpopulate the Fermi limits [11] and such a 
detailed study with different model parameters would have important multi-messenger im¬ 
plications. However, such a limit would also involve more parameters as one would need to 
understand the cascading of high energy gamma rays to the low energy flux and hence more 
uncertainties in the already parameter abundant model. The constraints coming only from 
the future IC neutrino data would be independent of uncertainties connecting the diffuse 
neutrino background with the diffuse gamma ray flux and should give the cleanest limits. 

The present analysis shows the importance of the multi energy study involving the SNRs. 
We focus on the HNRs as the higher energy counterpart. However, there are several other 
possibilities from relativistic CRB remnants to slow jet AGNs. The relative population and 
SFR of these rare, high energy cosmic accelerators would need more understanding. On the 
other hand the nature and the rate of the normal SNRs are better understood and the diffuse 
neutrino flux from SNRs is more natural compared to the extreme high energy accelerators 
contribution. Hence the SNRs contribution up to hundreds of TeV energies should give the 
reference and baseline for the fluxes of all the connected higher energetic astrophysical events. 
Note added: A week after the submission of the present work, [41] proposed a similar idea 
to the one presented here. 
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